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Abstract. The paper addresses generalized Borel summability of "]+" dif- 
ference equations in "critical time" . We show that the Borel transform Y of a 
prototypical such equation is analytic and exponentially bounded for Sft(p) < 1 
but there is no analytic continuation from toward +oo: the vertical line 
i := {p : 5R(p) = 1} is a singularity barrier of Y. 

There is a unique natural continuation through the barrier, based on the 
Borel equation dual to the difference equation, and the functions thus obtained 
are analytic and decaying on the other side of the barrier. In this sense, the 
Borel transforms are analytic and well behaved in C \ I. 

The continuation provided allows for generalized Borel summation of the 
formal solutions. It differs from standard "pseudocontinuation" UJ. This 
stresses the importance of the notion of cohesivity, a comprehensive exten- 
sion of analyticity introduced and thoroughly analyzed by Ecalle. 

We also discuss how, in some cases, Ecalle acceleration can provide a pro- 
cedure of natural continuation beyond a singularity barrier. 



1. Introduction 

In the case of generic differential equations, generalized Borel summation of a 
formal power series solution, in the sense of Ecalle U, essentially consists in the 
following steps: (1) Borel transform with respect to a critical time, related to the 
order of exponential growth of possible solutions, (see also the note below), usual 
summation of the obtained series, analytic continuation along the real line or in its 
neighborhood, proper averaging of the analytic continuations (e.g. medianization) 
toward infinity, possible use of acceleration operators and Laplace transform C. 

The choice of the critical time, or of a very slight perturbation -weak acceleration- 
of it is crucial for Ecalle summability. A slower variable (time) would hide the 
resurgent structure encapsulating the Stokes phenomena, and, perhaps more im- 
portantly, introduces superexponcntial growth preventing Laplace transformability 
at least in some directions. In a faster variable, convergence of the Borel trans- 
formed series would not hold. 

In some functional equations and so called type 1 + difference equations, new 
difficulties occur. For them, Ecalle replaces analyticity with cohesivity 5 . This 
property was studied rigorously for some classes of difference equations by Immink 
6 . It is the purpose of this note to show the importance of this notion: even in 
simple 1 + difference equations it is shown that critical time Borel transform has 
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barriers of singularities, preventing continuation in some half-plane. This occurs in 
the prototypical equation 



(example 2. of [HI ) - A simple proof of Borel space natural boundaries is not present 
in the literature, as far as the author is aware. We also show that the barrier 
is traversable: on the real line the associated function is well defined and Laplace 
transformable to a solution of the difference equation. This function is real analytic 
except at one point and, in fact has analytic continuation in the whole of C\£ with 
£ = {p : = 1} a singularity barrier. The present approach is adaptable to more 
general equations. 

We expect barriers of singularities to occur quite generally in 1 + cases, due to 
the fact that the pole position is periodic in the original variable, while critical time 
introduces a logarithmic shift in this periodicity. This leads to lacunary series in 
Borel plane, hence to singularity barriers. 

Nonetheless, further analysis shows that, in this simple case, and likely in quite 
some generality, softer Borel summation methods and study of Stokes phenomena 
are possible, relying on the convolution equation for continuation through singular- 
ity barriers. 

In spite of its simplicity, the properties in Borel plane of this equation, in the 
critical time, are very rich. 

Note on critical time. The solution of the homogeneous equation associated to 
JU, f(x) — l/T(x) has large x behavior (x/2Tr) 1 ^ 2 e~ xlnx+x . The critical time z is 
then the leading asymptotic term in the exponent, z = xhix [§]. (The origin of 
the terminology 1 + is related to the exponential order slightly larger than one of 
/). Various slight perturbations of this variable, weak accelerations, are used and 
indeed are quite useful. 



Theorem 1. Let Y(p) be the Borel transform of y in £7J) in the critical time z. 
Then Y(p) is analytic on {p ^ : arg(p) G — 2tt,tt + 27r);3i(p) < 1} and 
exponentially bounded as \p\ — > oo in this region. The line £ = {p : 5Rp = 1} is a 
singularity barrier ofY. 

Proof of the theorem. Let y be the formal power series solution of JJJ . We study 
the analytic properties of the Borel transform By := Y(p) of the on So, the Riemann 
surface of the log at zero, with respect to the critical time z. In critical time the 
functional equation of By © is unwieldy, and instead we look at the meromorphic 
structure of solutions on which we perform a Mittag-Lefner decomposition. 

It is straightforward to check that y is the asymptotic series for arg(x) ^ of 
the following actual solution of 



(1) 



y{x + 1) = -y(x) + - 



2. The singularity barrier 




oo 



(2) 



fc=l j=l 
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The fact that Hes(yo;x = n) = e~~ 1 /T(n) and the behavior at infinity of yo show 
that the Mittag-Leffler partial fraction decomposition of © is 

(1) Analyticity in the left half plane. The inverse function z i— » x(z) of xhix 
is analytic on So \ (— e _1 , 0) as it can be seen from the differential equation ^| = 
(1 + lnx) -1 . Then Y(p) is the analytic continuation of the function defined for p 
negative by 

(4) / ef a y (x(z))dz = ^ f ^"yo{x(z))dz, peW 

where C is a contour from oo + iO around — e _1 and to oo — iO. 

(2) Identities for finding continuation in {z : 3t(z) < 1} and exponential bounds. 
For analytic continuation clockwise we start from axgp = 7r and rotate up the 
contour, collecting the residues: 



1 V- 1 f ePZdz r^r n 1 /" V- 1 



2e7ri j^— ' T(fc) J c — fc 2e7ri j£~J r(fc)(ir(z) — fc) 



(5) where := ^ 



■° 1 + lnfc 



gpfc In /c 



and where for small > 0, Ci is the contour from ooe , ^ >+l ° around (— e _1 ,0) 
to ooe l ^~*°. As argp is decreased from to zero (and further to — 7r), can be 
increased from + to 2tt~ making J c visibly analytic in {p ^ : argp £ ( — 7T, 7r)} 
and exponentially bounded as \p\ — > oo. We decomposed Y into a sum of a lacunary 
Dirichlet series and a function analytic in the right half plane. 

(2) The natural boundary. The Dirichlet series F is manifestly analytic for -ftp < 
1. As p | 1 we have F(p) — > +oo and thus F is not entire. But then, by the Fabry- 
Wennberg-Szasz-Carlson-Landau theorem |S] pp. 18, I is a singularity barrier of F 
and thus of Y. For a detailed analysis, see also the note below. □ 



Note: Description of the behavior of F at I. Since all terms of the Dirichlet 
series are positive on the real line, it is easy to check using discrete Laplace method 1 
that F increases like an iterated exponential along K + toward £, F(p) oc exp((l — 
p) exp(l/(l— p))). There are densely many points near I where the growth is similar; 
it suffices to take a sequence of k £ N, 5ft (p) = k/(l + ln(fc)) and (1 + ln(fc))3(p) 
very close to an integer multiple of 2ir. (A Rouche type argument shows there are 
also infinitely many zeros with a mean separation of order the reciprocal of the 
maximal order of growth, ln(d) ~ —(1 — p)e 1 ^ 1 ~ p \) Rather than attempting some 
form of continuation through points where F is bounded, which are easy to exhibit, 
we prefer to soften the barrier first, by acceleration techniques. 



Determining, for fixed p, the maximal term of the series and doing stationary point expansion 
nearby. 
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3. General Borel summability in the direction of the barrier. 
Properties beyond the barrier. 

Strategy of the approach. It is convenient to perform a "very weak acceleration" 
to smoothen the behavior of Y(p) near i. The natural choice of variable is z — 
lnr(fc), but we prefer to slightly accelerate further, to z m (x) defined in Remark ^ 
below. We construct actual solutions of 0) starting from an incomplete Borel sum. 
We identify these actual solutions and show they are inverse Laplace transformable. 
Furthermore, they solve the associated convolution equation in Borel space. From 
these points of view, we have a unique continuation on R + . We show that the 
function thus obtained is real analytic on R \ {1} and continuable to the whole of 

The general solution of Q is 

(6) y(x) = y (x) + ^| 

where / is any periodic function of period one, as it can be easily seen by making 
a substitution of the form |JBJ in the equation. It can be easily checked that the 
following solution of Q 

ir cot irx 

(7) yi{x)=yo + - — 

is an entire function, and has the asymptotic behavior y, the formal series solution 
to defined in the proof of the theorem. 

Remark 1. Let m S N and z m (%) = x Inx— x—(m+^) hix. For given C > 0, there 
is a one-parameter family of solutions of 0) which are analytic and polynomially 
bounded in a region of the form Sc = {x : 5R(z m (x)) > C}. They are of the form 
y c (x) — yi(x) + c/r(x) for some constant c. 

Proof. The solution Q already has the stated boundedness and analyticity proper- 
ties (and in fact, it decreases at least like x~ m in Sc)- The general solution is of the 
form y\ + f(x)/T(x) with / periodic, as remarked at the beginning of the section. 
Analyticity implies / is analytic and boundedness in the given region implies / is 
bounded on the line dSc- By periodicity, / is polynomially bounded in the whole 
of C, which means / is a polynomial, and by periodicity, a constant. □ 

Theorem 2 (Generalized Borel summability). (i) There exists a one parameter 
family of solutions of 0) which can be written as L Zm H c :— J Q e~ ZmP H c (p)dp 
where H c = B Zm y is analytic and exponentially bounded for 3?(p) < 1 and H c g 
C m - 1 (M + ). 

(ii) H c are real analytic on R + \ {1}; they extend analytically to <C\£, and £ is a 
singularity barrier H c and the functions are C m_1 on the two sides of the barrier 2 . 
Furthermore, for 3?(p) > 1, H c decrease toward infinity in C. 

Remark 2. It would not be correct at this time to conclude that, say, C~ 1 y\ pro- 
vides Borel summation of y; we need to show that y\ satisfies the necessary Gevrey- 
type estimates to identify the inverse Laplace transform with By in the unit disk. We 
prefer to proceed in a more general way, not using explicit formulas, but construct- 
ing actual solutions starting with an incomplete Borel summation (and identifying 
them later with the explicit formulas). 

2 The values on the two sides cannot, obviously, be the same. 
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Proof of Theorem^ (i) We redo the analysis of the proof of Theorem 1 in the 
variable z — z m and we get a decomposition of the form J^J), where now F is 
replaced by 

00 In h A- BL r , 1 

^ p 2= y^ ^ + fc c p[fclnfc-fc-(m+^)lnfcJ 

k=l eF(fc) 

which is a Dirichlet series of the same type as F and hence has £ as a singularity 
barrier. However, F2 is (manifestly) uniformly C"™^ 1 up to £ and so is thus Y{p). 

For the solutions of that decrease in a sector in the right half -plane it is 
clear that the dominant balance is between y{x + 1) and 1/x. We then rewrite the 
equation to prepare it for a contraction mapping argument in Borel space. By a 
slight abuse of notation we write y(z) for y(x(z)) and we have 

(x{z) - l)y(x(z)) = y{x(z) - 1) + 1 

(x(z) - l)y(z) = y(z - g(z)) + 1 
where g(z) = lnz — lnlnz + o(l) and then 

(x{z)-l)y{z) = Y,y^{z)g(z) k /k\ + l 

k=0 

Thus, dividing by x(z) — 1 and taking inverse Laplace transform, with Gk(p) the 
inverse Laplace transform of g(z) k / (x(z) — l)/k\, we have 

00 

(9) Y(p) = J2l(-p) k Y}*G k (p)+F(p) 

fc=0 

The term Gk is (roughly) bounded by |e~ fe ( 1-p )|, as can be seen by the saddle 
point method applied to the inverse Laplace transform integral. It is easy to check, 
using standard contraction mapping arguments (see e.g. |2]), that Y is given by a 
convergent ramified expansion in the open unit disk. This was to be expected from 
estimates of the divergence type of the formal solutions of However, given the 
estimates on the terms of the convolution equation, the equation, as written, cannot 
be straightforwardly interpreted beyond 5i(p) = 1, the threshold of convergence of 
the ingredient series. It is however possible to write a meaningful global equation 
by returning to the definition in terms of Laplace transform. We then write 

1 />c+ioo poo poo 

C- l y{z + g{z)) = — I dze pz dqe- q{z+9[z)) Y{q) = H{p,q)Y(q)dq 

27H Jc^ioc Jo Jq 

where 

-1 />e+zoo -1 />c+zoo 

H(p,q) = : / e (P-l)z-qg(z) dz = / e (p-g)z+?(lnln 2 + ...) z -g rfz 

2tT* Jc—ioo 27TI J c _,; 00 

which is well defined for q > and integrable at q = 0; the convolution equation 
becomes 



/•OO 

(10) / H(p,q)Y(q)dq = Y*£- 

Jo 



r 1 


+ C- 1 


f 1 ] 


-x(z) — 1- 


-x(z) — 1- 





Based on the solution on [0, 1) of © we construct solutions to and their 
inverse Laplace transforms provide continuation of Y past 5R(p) = 1 and implicitly 
solutions to (|TU|l . 
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Wc define the incomplete Borel sum 




e~ zp Y 1 (p)dp 



Formal manipulation shows that y satisfies with errors of the form 3 o(e~ z ) or 
o(x m /T(x)) in the variable x) where the estimate of the errors is uniform in the 
right half-plane in z, or in a region Sc w.r. to x. 

We look for a solution of in the form y + S(x)/T(x). Then S(x) satisfies S(x + 
1) = S(x) + R(x) (the 1 + degeneracy is not present anymore) where R(x) = o{x m ) 
with differentiable asymptotics (by Watson's lemma). A solution of this equation is 
-P m+3 Y^kLo -R*- m+3 ' ) (x + k) , with V an antiderivative, which is manifestly analytic 
and polynomially bounded in regions of the form Sc, and y + S/T is manifestly a 
solution of P|l. which, by construction, is also polynomially bounded in Sc- 

By Remark H V + S/T is one of the solutions y c . But y c is inverse Laplace 
transformable with respect to z, and has sufficient decay to ensure the existence of 
m — 1 derivatives of the transform. By Remark ^ any solution that decreases in 
the natural region Sc in the right half plane can be represented in this way and 
thus the conclusion follows. □ 

Corollary 3. In {p : 5i(p) < 1}U [1, oo), there is a one parameter family of Laplace 
transformable solutions to fl(J\). the functions H c in Theorem^ (i). They have I as 
a barrier of singularities. 

Proof of Theorem^ (ii). Since all Laplace transformable solutions to l|l(J|) are those 
provided in Remark ^ we analyze the properties of the inverse Laplace transform 
of these functions for 3?(p) > 1. 

We note that, due to the fact that y c (z m ) increase at most as e Zm /z™, we can 
deform for 3?(p) > 1, the integral 



where C starts at — oo — ie, avoids the origin through the right half plane and turns 
back to — co + ie. In view of the bound mentioned above for y c (z m ), this function is 
manifestly bounded and analytic for 3?(p) > 1, and in fact is continuous with m— 1 
derivatives up to 5R(p) = 1. 

Cohesive continuation and pseudocontinuation. It follows from our analysis 
and from the fact that Ecalle's cohesive continuation also provides solutions to the 
equation, that the results of the continuations are the same (modulo the choice of 
one parameter, discussed in the Appendix). This type of continuation is the natural 
one since it provides solutions to the associated convolution equation. It is easy 
to see however that this continuation is not a classical pseudocontinuation through 
the barrier, as it follows from the following Proposition. 

Proposition 4. The values of H c on the two sides of£ are not pseudocontinuations 
[HI of each-other. 



(11) 




771 



to an integral 



(12) 




^Resulting from incomplete representation of l/(x(z) — 1). 
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Proof. Indeed, pseudo-continuation [Hj, pp. 49 requires that the analytic elements 
coincide almost everywhere on the two sides of the barrier. But H c is continuous on 
both sides, and then the values would coincide everywhere, immediately implying 
analyticity through £, a contradiction. 

Remark 3. The axis M + , which is also a Stokes line, plays a special role. No other 
points on the singularity barrier can be used for Borel summation, as shown in the 
proposition below. 

Proposition 5. No Laplace transformable solution of HKJp exists, in directions 
e z ^M. + , 4> s (0, 7r/2). (The same conclusion holds with <fi S (— 7r/2, 0).) 

Proof. Indeed, the Laplace transform y of such a solution would be analytic and 
decreasing in a half plane bisected by e 1 ^ and solve©. Since l/r(a;) is entire and 
the general solution is of the form © , by periodicity /i = / — - cot irx would be 
entire too. Taking now a ray te 1 ^^ 7 '/ 2 "^ we see, using again periodicity, that f\ 
decreases factorially in the upper half plane. Standard contour deformation shows 
that half of the Fourier coefficients are zero, f\(x) — J2keN c ke tkx and that, because 
/ is entire, c& decrease faster than geometrically. But then fi(x) =: F (exp(2mx)) 
with 1 1 ► F(t) entire. When x — > ioo, t — > and, unless F — 0, we have F(t) ~ ct n 
for some n € N, thus f(x) ~ ce mx , incompatible with factorial decay. This means 
/ = but then © is not analytic on the real line 4 . □ 

4. Appendix: Weak acceleration, integral representation, median 
choice, natural crossing of the barrier 

A weak acceleration is provided by the passage xlnx — x i— > x. The x- inverse 
Laplace transform of satisfies e~ p Y ~ f£ Y(s)ds — 1 = with the solution 
Y = e _1 exp(p + exp(p)). CY exists along any (combination of) paths R n starting 
from the origin and ending on a ray of the form p = K + + (2n + n G Z. The 
function / + = J R e~ xp e p+e ~ x dp is manifestly entire 5 . For x = —t\ t — > oo the 
saddle point method gives 

/ — v /2^ e *lnt-t+7rit+i lnt-l 

which identifies / + with y\ + m/e/T(x). With obvious notations, we see that 
yi = + /-), reminiscing of medianization. We have also checked numerically 
that yi is approximated by least term truncation of its asymptotic series with errors 
o(l/r(x)). (The integral representation would allow for a rigorous check, but we 
have not done this and we state the property as a conjecture; we also conjecture that 
the solution constructed in Proposition is yi; this could be checked by looking 
at the asymptotic behavior on dSc-) There is, obviously, only one solution so 
well approximated. It should then be considered as the natural candidate for the 
medianized transform in critical time and its inverse Laplace transform, defined on 
the whole of M + , and the natural continuation of the Borel transform By past the 
barrier. For all these reasons it is likely, but we have not checked it rigorously, that 
yi corresponds to the medianized cohesive continuation of Ecalle. 

4 We should note that a procedure mimicking the proof of Theorem |5] (i) in non-horizontal 
directions would fail because now the remainders R(x) would grow fast along the direction of 
evolution — parallel to R+ . 

5 It provides, in view of the superexponential properties of the integrand, Borel oversummation. 
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Remark 4. The procedure described of naturally crossing a barrier does not neces- 
sarily depend on the existence of an underlying functional equation. It is sufficient 
to have accelerations as above that allow for Borel (over) summation along some 
paths, and choose as a natural actual function the one that has minimal errors 
in least term truncation or resort to a medianized choice. The process of contin- 
uation through the barrier can be written as the composition C~^L zl B zl L Zm with 
L formal Laplace transform, and is expected to commute with most operations of 
natural origin. It is applicable to many other series including the Dirichlet series 

V°° pfP- 1 )" 2 
Z^fe=0 e 

Finally, it seems a plausible conjecture that in the case of nonlinear systems, 
infinitely many equally spaced "isolated" barriers should occur. 
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